V 


C  -- 


i 


'-1 


Williams/Garrison/Mourad 


2  April  1992 

Experimental  examination  of  growing  and  newly  submerged  sea  ice  including 
acoustic  probing  of  the  skeletal  layer 


K.L.  Williams,  G.R.  Garrison,  and  P.D.  Mourad 

Applied  Physics  Laboratory,  College  of  Ocean  and  Fishery  Sciences,  University  of  Washing¬ 
ton,  Seattle,  Washington  98105 

(Received 

Results  of  an  in-situ  experiment  to  examine  acoustic  parameters  and  scattering 
properties  of  growing  and  newly  submerged  arctic  sea  ice  are  presented.  The  pri¬ 
mary  emphasis  is  on  the  acoustic  properties  of  highly  porous  sea  ice.  High  resolu¬ 
tion  (to  about  1  cm)  vertical  profiles  of  the  longitudinal  wave  speed  in  growing  and 
submerged  sea  ice  are  derived  from  the  experimental  data.  The  results  indicate  a 
skeletal  layer  about  3  cm  thick  at  the  bottom  of  the  growing  ice.  The  time  depen¬ 
dence  of  venical  sound  speed  profiles  in  a  submerged  ice  block  reveals  a  long  time 


OTIC 

^  ELECTE 

^  1992 

CNJS 

mg 

mS 

CM^ 


I 

D 

< 


scale  process  of  up  to  80  hours  and  a  short  process  of  a  few  hours.  The  former  pro- 

t- 

cess  is  the  warming  of  the  block;  the  latter  is  thought  to  be  due  to  the  disturbance  of 
hydrostatic  equilibrium  as  the  ice  is  submerged.  Concurrently  with  the  wave  speed 
measurements,  scattering  from  the  undersurface  of  the  ice  was  measured  at  several 
frequencies.  A  comparison  with  predictions  based  on  the  sound  speed  data  demon¬ 
strates  the  ability  to  predict  normal  incidence  ice  reflectivity  from  sound  speed 
profiles  as  well  as  the  viability  of  using  scattering  data  in  inversions  to  obtain  sound 
speed  profiles.  Ab.sorption  of  a  longitudinal  wave  propagating  venically  in  ice  was 
also  measured.  The  peak  absorption  rate  found  in  the  skeletal  layer  was  between  2 
and  5  dB/cm  at  92  kHz.  The  temperature  dependence  of  absorption  seen  in  the  sub¬ 
merged  ice  suggests  that  the  McCammon-McDaniel  equation  (~T  is  useful 
away  from  the  skeletal  layer.  The  measurements  of  ice  temperature  as  a  function  of 
time  and  depth  allow  a  calculation  of  the  themtal  diffusiviiy  of  the  upper  region  of 
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growing  sea  ice  (0.0080  cm^/s)  and  an  “effective”  thermal  diffusivity  of  the  sub¬ 
merged  ice  block  (about  0.001 5  cm^/s).  Measurements  of  salinity  profiles  in  the  ice, 
along  with  the  temperature  profiles,  are  used  to  calculate  porosity,  which  is  then 
used  in  models  of  the  longitudinal  wave  speed  as  a  function  of  porosity.  A  compar¬ 
ison  between  models  and  experiment  is  employed  in  suggesting  a  possible  mecha¬ 
nism  for  the  increased  speed  seen  experimentally  in  the  bottom  portion  of  the  ice 
after  submergence. 

PACS  numbers:  43.20.Fn,  43.30.Ma,  43.30.Pa 
INTRODUCTION 

The  acoustic  properties  of  sea  ice  are  inherently  important  to  studies  of  sound  propagation  and 
scattering  in  the  Arctic.'  '  The  richness  of  the  acoustic  behavior  of  sea  ice  stems  from  its  saline 
nature.  In  panicular,  sea  ice  is  the  result  of  the  solidification  of  a  sub-eutectic  melt. This  accounts 
for  the  high  porosity  layer  at  the  ice/water  interface,  known  as  the  “skeletal”  layer.  The  high  po¬ 
rosity  makes  measurements  of  the  acoustic  properties  in  the  skeletal  layer  extremely  sensitive  to 
penurbations.  Because  this  layer  is  a  few  centimeters  thick,  it  is  of  primary  importance  for  model¬ 
ing  the  scattering  of  high  frequency  (above  10  kHz)  acoustic  energy  from  the  underside  of  the  ice. 

The  porosity  of  sea  ice  is  controlled  by  its  salinity  and  temperature.  ^  For  a  constant  salinity, 
an  increase  in  temperature  is  accompanied  by  an  increase  in  porosity.  Thus  if  a  block  of  sea  ice  is 
submerged,  and  thereby  warmed  during  keel  formation,  the  entire  block  may  become  highly  po¬ 
rous.  In  this  case  the  dimensions  of  the  block,  or  of  a  structure  made  from  many  blocks,  are  the 
appropriate  length  scales  to  use  in  determining  the  frequency  range  where  the  acoustic  behavior  of 
high  porosity  ice  is  important.  Since  the  ice  blocks  that  fomi  keels  in  the  Arctic  have  length  scales 
on  the  order  of  meters  and  the  keels  h:ive  .scales  on  the  order  of  tens  of  meters,  that  frequency  range 
extends  below  100  Hz. 

Because  high  porosity  sea  ice  is  sensitive  to  disturbance,  we  measured  its  acoustic  properties 
in  situ.  We  also  measured  salinity  and  temperature,  not  only  for  information  about  the  physical  pa¬ 
rameters  of  the  ice  (e.g.,  porosity  and  thennal  diffusivity),  but  to  relate  these  properties  to  the  mea- 
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sured  acoustical  properties  (longitudinal  wave  sound  speed,  reflectivity,  and  absorption).  The 
experimental  apparatus  and  method  have  produced  high  vertical  resolution  measurements  (to  about 
1  cm)  of  acoustic  properties  in  high  porosity  sea  ice  both  during  growth  and  after  submergence. 

The  anicle  is  organized  as  follows.  Section  I  describes  the  experimental  instrumentation,  ar¬ 
rangement,  and  procedure,  as  well  as  some  of  the  environmental  conditions  during  the  experiment. 
Section  II  gives  the  main  results,  dealing  first  with  growing  ice  and  then  with  submerged  ice.  Sec¬ 
tion  III  is  a  summary  of  results.  The  appendix  details  the  calculations  for  obtaining  the  thermal  dif- 
fusivities  quoted  in  the  body  of  the  paper. 


I.  EXPERIMENT 

The  experiment  was  conducted  at  an  Applied  Physics  Laboratory  Ice  Station  (APLIS)  estab¬ 
lished  in  the  Beaufort  Sea  and  occupied  from  the  end  of  February  to  the  middle  of  April  1990.  In 
the  first  part  of  the  experiment,  a  2-m  square  hole  was  cored  in  the  ice  and  sensor  arrays  were  sus¬ 
pended  in  the  water  to  monitor  conditions  as  ice  formed  around  them.The  ice  grew  to  about  42  cm 
thickness  in  10  days.  During  this  time,  the  temperature  and  salinity  of  the  growing  ice  were  mon¬ 
itored  while  measurements  were  made  of  the  vertical  longitudinal  wave  speed  and  attenuation.  The 
normal  incidence  acoustic  reflectivity  of  the  growing  sea  ice  was  also  measured  from  a  transducer 
below  at  frequencies  from  40  to  220  kHz. 

In  the  next  phase  of  the  experiment,  an  ice  block  was  cored  from  a  ponion  of  the  newly  grown 
ice  containing  the  sensors.  This  block,  a  cylinder  about  58  cm  in  diameter  and  42  cm  in  height,  was 
submerged  and  held  in  place  while  we  made  periodic  measurements  of  temperature,  longitudinal 
wave  sound  speed,  attenuation,  and  reflectivity  (similar  to  those  taken  in  the  first  pan  of  the  exper¬ 
iment)  until  the  ice  reached  thermal  equilibrium  with  the  surrounding  water  (about  80  hours). 

Details  of  the  experiment  are  given  in  the  following  subsections.  Section  l.A  discusses  the  sen-  „ 

sor  arrays  and  the  transmitter/receiver  acoustic  transducers.  Section  l.B  outlines  the  experiment  ar- 
rangeinent  and  procedure.  Section  I.C  describes  some  general  environmental  conditions.  ;  ** 
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A.  Instrumentation 


The  in-water  instrumentation  comprised  two  vertical  thermistor  arrays,  a  venical  hydrophone 
array,  and  two  transmitter/receiver  acoustic  transducers.  The  on-ice  instrumentation  consisted  of 


switching  boxes  for  selecting  particular  sensors  in  the  vertical  arrays.  All  electronics  control,  signal 
generating,  and  data  recording  instrumentation  were  housed  in  a  plywood  building  about  50  m 
from  the  experiment  site. 

The  two  transmitter/receiver  units  were  a  “platter”  and  a  “sonar”  transducer.  The  platter^  was 
used  for  the  normal  incidence  reflectivity  measurements  and  to  measure  ice  growth.  An  indepen¬ 
dent  measure  of  ice  growth  was  made  with  the  sonar^  transducer. 


The  thermistor  and  hydrophone  arrays  were  built  in-house.  Two  thermistor  arrays  were  con¬ 
structed  using  Yellow  Springs  Instrument  (YSI)  preci.sion  thermistors  with  impedance-to-temper- 
ature  calibrations  matched  to  0.2°C.  Each  thermistor  was  encased  in  marine  epoxy.  The 


hydrophone  array  was  made  of  PZT  ceramic  elements  approximately  cubic  in  shape  (0.5  cm  per 
side)  with  a  resonance  frequency  around  200  kHz.  Each  element  was  placed  at  the  end  of  a  small 
Plexiglas  tube  (0.5  cm  o.d.  and  5  cm  in  length)  and  encase"^  in  PRC  (Products  Research  and  Chem¬ 
ical  Company)  potting  compound.  The  arrays  were  constructed  using  Plexiglas  tubing  (o.d.  1.5-2 
cm  and  length  about  90  cm)  and  marine  epoxy.  The  wire  leads  running  from  arrays  to  switching 
boxes  were  enclosed  in  surgical  tubing.  Each  thermistor  array  contained  24  sensors  spaced  approx¬ 
imately  2  cm  apart.The  hydrophone  array  had  12  sensors,  the  upper  half  spaced  about  5  cm  apart 
and  the  lower  half  about  2.5  cm  apart.  The  total  vertical  coverage  of  the  arrays  (45  cm)  was  based 
on  predicted  freezing  rates  and  the  time  allotted  for  the  experiment  (two  weeks). 

Figure  1  shows  two  of  the  vertical  arrays  (hydrophones  to  the  left).  Note  that  the  hydrophones 
were  arranged  in  a  helix  to  prevent  acoustical  shadowing.  Figure  2  gives  the  vertical  spacings  of 
these  sensors  and  the  horizontal  spacing  between  the  thennistor  and  hydrophone  arrays. 

Switching  boxes  on  the  ice  were  connected  to  Nim-^in  control  modules  built  in-house.  The  hy-  Hith 
drophone  control  module  allow-ed  access  to  any  2  of  the  12  hydrophones.  The  thermistor  control 
module  accessed  any  .sensor  in  the  thermistor  arrays.  Wavefomi  genertition,  timplification.  and  re- 
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FIG  1.  Array  of  hydrophones  (at  left)  and  thermistors  (at  right)  that  were  suspended  in  the  center 
of  a  2-m  square  hole  in  the  ice.  Both  arrays  were  suspended  from  a  wood  2x4  across  the  hole. 
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FIG  2.  Depths  of  the  suspended  hydrophones  and  thermistors.  The  top  thermistor,  #24,  was  at  the 
ice/air  interface.  Inoperable  sensors  (e.g.,  #13)  have  been  omitted. 
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ception  for  the  platter  transducer  were  performed  by  the  APL  MAME^®  system.  The  sonar  trans¬ 
ducer  had  its  own  control  box.  All  acoustic  signals  were  recorded  on  a  Nicolet  digitizing  scope. 
The  resistances  measured  for  each  thermistor  were  recorded  by  hand. 

B.  Experimental  Arrangement  and  Procedure 

The  ice  canopy  was  about  2  m  thick  at  the  experiment  site.  Two  2-m  square  holes  were  cut 

through  the  ice  canopy,  one  for  the  main  experiment  and  the  other  for  monitoring.  The  later  was  3 
m,  center  to  center,  from  the  main  hole.  One  of  the  thermistor  arrays  and  the  hydrophone  array  were 
vertically  supported  at  the  center  of  the  main  hole.  The  other  thermistor  array  was  suspended  at  the 
center  of  the  monitor  hole.  Each  thermistor  array  was  arranged  so  that  the  top  sensor  was  at  the 
water/air  interface.  Figure  3  is  a  diagram  of  the  in-water  instrumentation  in  the  main  hole  during 
the  submergence  phase,  discussed  later. 

Salinity  measurements  were  taken  on  sections  of  ice  cores  from  the  monitor  hole  during  the  ice 
growth  phase.  Cores  were  taken  from  this  hole  to  avoid  disturbing  the  main  hole  from  which  an  ice 
block  was  to  be  cored  for  the  second  part  of  the  experimeni.The  second  thermistor  array  let  us  de¬ 
termine  whether  ice  samples  in  the  two  holes  were  in  equivalent  states.  The  effort  to  obtain  cores 
was  worthwhile,  considering  the  importance  of  salinity  in  determining  ice  properties  (e.g.,  poros¬ 
ity).  However,  the  invasive  nature  of  the  technique  causes  the  salinity  measurements  to  be  of  lower 
resolution  and  quality  than  the  temperature  and  acoustic  parameter  measurements. 

The  platter  transducer  shown  in  Fig.  3  was  suspended  from  three  lines  120°  apart,  which  were 
easily  adjusted  to  place  the  transducer  directly  below  the  hydrophone  array.  The  “sonar"  trans¬ 
ducer  was  on  a  rotating  arm  so  that  it  could  be  placed  below  the  arrays  during  the  ice  growth  phase 
and  then  moved  away  before  the  ice  block  was  submerged.  The  short  acoustic  pulses  emitted  from 
these  transducers  were  reflected  from  the  ice/  water  interface  and  then  received  back  at  the  trans¬ 
ducers.  They  were  used  to  determine  ice  thickness  as  a  function  of  time,  as  well  as  the  reflectivity 
of  the  ice  as  a  function  of  frequency.  The  platter  transducer’s  emissions  were  received  by  the  h\- 
drophone  array,  allowing  measurement  of  ice  acoustic  parameters 

Measurements  for  determining  reflectivities,  acoustic  parameierN,  and  temperature  of  the  ice 
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FIG.  3.  Arrangement  after  the  ice  in  the  2-m  hole  had  grown  to  42  cm  and  a  block  of  this  new  ice 
had  been  cored  and  submerged.  (Depth  of  transducer  not  to  scale.) 
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were  made  several  times  each  day  during  the  ice  growth  phase  of  the  experiment.  After  submer¬ 
gence  of  the  ice  block,  the  same  type  of  measurements  were  made,  but  more  frequently  (30  mea¬ 
surements  in  the  80  hours  monitored  after  submergence).  The  results  derived  from  these 
measurements  and  the  salinity  measurements  form  the  bulk  of  Sec.  II. 

C.  Environmental  Conditions 

A  portable  weather  station  was  operated  at  APLIS  to  record  air  temperature,  wind  speed,  etc. 
Figure  4  shows  a  curve  of  the  air  temperature  recorded,  along  with  the  temperatures  determined 
from  the  top  sensors  in  the  two  thermistor  arrays  during  the  ice  growth  phase.  The  array  tempera¬ 
tures  tended  to  be  higher  than  the  weather  station  temperatures.  This  could  be  due  to  snow  coverage 
when  winds  cause  drifting.  The  holes  were  cleared  of  snow  periodically,  but  not  always  immedi¬ 
ately  before  measurement.  Note,  however,  that  the  two  thermistors  are  always  in  close  agreement, 
which  indicates  that  the  ice  freezing  process  was  being  driven  by  the  same  thermal  gradient  at  the 
two  holes.  Note  too  that  in  the  last  day  of  the  ice  growth  phase  a  significant  front  moved  in,  with 
air  temperature  increasing  more  than  10°.  This  rapid  rise  allowed  a  calculation  of  thermal  diffusiv- 
ity  using  the  vertical  profiles  of  the  temperature  of  the  ice  during  that  time. 

In  Fig,  5  the  salinity  and  temperature  profiles  measured  from  ice  cores  extracted  from  the  42- 
cm  ice  in  the  monitor  hole  on  the  last  day  of  the  ice  growth  phase  are  compared  with  profiles  for  a 
region  of  the  ice  canopy  near  the  experiment.  Data  for  the  deepest  few  centimeters  of  high  porosity 
ice  are  not  included  because  brine  drainage  during  extraction  typically  makes  them  unreliable.  The 
computed  porosity  and  the  longitudinal  wave  sound  speed  derived’^  using  this  porosity  are  also 
shown  in  this  figure.  The  low  salinity  region  near  the  top  of  the  canopy  profile,  and  the  overall  low 
salinity,  indicate  that  the  ice  in  this  region  has  gone  through  at  least  one  warming  period  and  is  thus 
probably  multiyear  ice.  In  Sec.  11.4,  speeds  predicted  from  Ref.  1 1  will  be  compared  with  those 
experimentally  determined  in  the  present  effon. 
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FIG.  4.  A  comparison  of  air  temperatures  from  the  camp  weather  station  with  temperatures  record¬ 
ed  by  thermistors  located  at  the  ice/air  interface.  Array  temperatures  tend  to  be  higher  than  weather 
station  temperatures,  probably  owing  to  snow'  cover. 
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FIG.  5.  Measured  temperature  and  salinity  of  ice  cores  taken  from  the  canopy  in  the  vicinity  of  the 
ice  growth  experiment  (dashed  lines).  The  porosit\-  and  sound  speed  profiles  are  calculated. The 
solid  lines  are  a  similar  presentatio*^  for  data  obtained  from  cores  taken  on  the  last  day  of  the  ice 
growth  experiment. 
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II.  RESULTS 
A.  During  Ice  Growth 

1.  Freezing  rate 

At  the  beginning  of  the  experiment,  with  the  platter  transducer  suspended  beneath  the  2-m  hole, 
we  measured  acoustic  reflections  from  the  surface  of  the  hole  after  the  water  surface  was  swept 
clear  of  ice  crystals.  Measurements  were  taken  at  several  frequencies:  37,  52, 92, 1  SO^and  220  kHz. 
Reflections  arriving  before  the  surface  reflection  were  identified  as  ei  her  the  bottom  of  the  ice  can¬ 
opy,  the  arm  that  supported  the  ..onar  transducer,  or  the  lower  end  of  the  sensor  strings.  These  fixed 
references  were  necessary  because  the  platter’s  depth  varied  as  the  supporting  lines  stretched. 

For  the  first  part  of  the  reflection  measurements,  ice  growth  was  determined  by  comparing  the 
reflection  time  for  the  ice  with  the  reflection  time  for  one  of  the  sensor  strings  at  220  kHz.  Later, 
the  rotating  arm  reflection  fell  within  the  digitizing  window  and  the  ice  reflection  time  was  then 
compared  with  the  arm  reflection  time  at  92  kHz.  In  determining  ’ce  growth,  we  used  a  sound  ve¬ 
locity  for  water  of  1436  m/s,  calculated  using  CTDs  taken  nearby.  Increase  in  ice  depth  for  the 
growth  period  is  shown  in  Fig.  6. 

Also  shown  are  the  results  of  an  independent  measure  of  ice  growth  using  the  sonar  transducer, 
which  operated  at  92  kHz  with  short  (0.5  ms)  cw  pulses.  The  ice  reflection  return  for  this  transducer 
was  recorded  on  a  digital  oscilloscope  about  twice  per  day.  To  obtain  the  ice  depth,  we  compared 
the  return  times  with  those  from  the  air/water  interface  at  the  start  of  the  e.xpe.'iment, 

A  curx'e  has  been  drtiwn  in  Fig.  6  to  best  represent  both  sets  of  data  It  gives  the  distance  from 
initial  water  surface  to  the  lower  surface  of  the  ice,  to  an  accuracy  of  about  0.5  cm. 

2.  Temperature  monitoring  and  calculation  of  thermal  diffusi\  ify 

Readings  of  the  thermistor  string  were  taken  three  to  six  times  per  day  during  the  ice  growtli 

phase  of  the  experiment.  The  measured  resistances  were  converted  to  temperature  with  a  calibra¬ 
tion  curve,  supplied  by  YSl  and  confimied  when  the  string  was  prepared  in  the  laboratory. 

A  plot  of  temperature  vs  time  for  each  thermisK'r  is  shown  in  Fig  7  The  thermistors  near  the 
surface  follow  the  air  temperature  (Fig.  4 ).  The  air  temi'ierature  c  ai  i.iiioiis  appe.ir  to  penetrate  about 
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FIG.  6.  Ice  depth  during  the  ice  growth  phase  as  determined  acoustically  from  two  different  trans¬ 
ducers  suspended  below  the  hole  Initially,  the  water  surface  was  the  reference  for  returns  to  the 
platter  transducer;  later,  because  the  platter  depth  varied,  the  reflections  from  the  bottom  of  the  sen¬ 
sor  array  and  the  am  for  the  sonar  transducer  were  used 
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FIG  7.  The  lemperatures  recorded  at  each  thermistor  during  ice  growth.  The  temperature  scale  has 
been  insened  from  that  for  the  air  temperatures  in  Fig  4  so  that  the  graphs  could  be  nested. 
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20  cm  into  the  ice,  to  thermistor  14,  with  a  lag  due  to  diffusivity.  Temperature  profiles  for  the  same 
data  are  plotted  in  Fig.  8.  The  first  profile  at  the  right  (No.  1)  shows  a  water  temperature  for  all 
thermistors  below  18  cm,  and  a  colder  temperature  for  all  above  1 8  cm,  indicating  that  the  ice  had 
frozen  to  that  depth  before  temperature  readings  began.  (There  was  a  delay  in  temperature  moni¬ 
toring  because  of  hardware  problems  in  the  thermistor  NIM-bin  module.)  With  time,  most  of  the 
temperatures  decreased  except  for  changes  related  to  air  temperature  variations.  On  29  March 
when  the  air  temperature  rose  to-14°C,  the  upper  30  cm  of  the  profiles  shifted  as  shown  by  the  four 
dotted  lines  in  Fig.  8.  The  coldest  profile  (No.  25)  was  for  0915  on  28  March,  the  time  of  the  coldest 
air  temperature  (Fig.  4). 

We  used  the  change  in  the  ice  temperature  profile  on  29  March  to  calculate  thermal  diffusivity 
assuming  a  constant  diffusivity  in  the  upper  region  of  the  ice  (the  ponion  below  -6°C).  The  calcu¬ 
lation  (detailed  in  the  appendix)  indicates  a  thermal  diffusivity  of  abot"  ^  0080  cm^/s.  In  reality, 
the  diffusivity  is  a  function  of  salinity  and  temperature,  with  the  most  severe  temperature  depen¬ 
dence  for  values  above  -6°  C.  This  becomes  evident  in  the  calculation  for  thermal  diffusivity  car¬ 
ried  out  for  the  submerged  ice  block  (Sec.  II. B).  Reference  12  gives  approximate  equations  for  the 
salinity  and  temperature  dependence  of  thermal  diffusivity  for  temperatures  above  -8°C.  These 
equations  predict  a  themial  diffusivity  of  0.(X)54  cnr/s  at  -8°C  (also  shown  in  the  appendix).  The 
upper  limit  of  the  diffusivity  can  be  set  using  the  nonsaline  ice  value  of  about  0,012  cm^/s  at  -20° 
C.  These  results  are  consistent  with  the  value  given  above. 

3.  Sound  speed  monitoring  and  initial  estimate  of  skeletal  zone  thickness 

Travel  times  between  hydrophones  were  measured  when  the  array  was  first  placed  in  the  water. 

These  times,  along  with  the  known  sound  speed  in  the  water  (C,^.),  give  the  separations  between 
the  acoustic  centers  of  the  hydrophones.  We  calculated  the  average  sound  speed  between  any  two 
hydrophones  at  various  times  after  they  were  frozen  in  by  dividing  the  spacing  by  the  travel  time. 
Although  we  did  not  have  enough  electronics  to  receive  each  ping  at  all  hydrophones,  we  were  able 
to  receive  at  two  adjacent  hydrophones  for  the  same  ping.  By  progressing  through  the  arra\ .  we 
obtained  data  th;ii  were  not  affected  by  ping-to-ping  variations  in  tlie  tr.insmitter  output. 
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F)G.  8.  Temperature  profiles  for  the  data  of  Fig  7.  Several  profiles  are  labeled  by  run  number;  the 
graph  at  the  top  of  the  figure  associates  run  number  u  ith  local  time.  The  dotted  lines  show  the  cffwCt 
of  a  12'C  rise  in  ambient  temperature  on  29  March 
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There  were  two  major  difficulties  with  the  experimental  data.  First,  hydrophone  10  had  a  very 
low  signal  to  noise  ratio  (SNR)  compared  with  the  other  hydrophones,  making  the  in-water  cali¬ 
bration  timing  measurement  difficult.  During  much  of  the  experiment  the  SNR  problem  was  com¬ 
pensated  for  by  using  the  Nicolet’s  averaging  capabilities.  The  second  difficulty  was  that  some  ice 
had  fomied  to  a  depth  below  hydrophone  1 1  before  the  “in-water”  measurements  of  travel  times 
were  made.  However,  good  estimates  of  the  separation  between  acoustic  centers  for  hydrophones 
10-12  were  obtained  by  comparing  all  calibration  times  with  rule-measured  spacings. 

The  sound  speed  in  each  hydrophone  interval  is  plotted  in  Fig.  9  as  a  function  of  time.  As  the 
ice  froze  in  each  interval  the  sound  speed  rose  above  After  the  speed  reached  about  3500  m/s, 
it  varied  (presumably)  because  of  changes  in  air  temperature  which  penetrated  the  ice  and  changed 
its  porosity. 

Figures  6  and  9  can  be  used  to  make  an  initial  estimate  of  the  thickness  of  the  “acoustic”  skel¬ 
etal  layer,  which  we  here  define  as  the  thickness  of  the  layer  between  the  ice/water  interface  and 
the  point  where  the  vertical  speed  reaches  3500  m/s.  We  chose  this  definition  because  variations 
in  speed  due  to  changes  in  atmospheric  conditions  are  apparent  at  higher  speeds  (Fig.  9). 

Consider  the  curve  in  Fig.  9  for  hydrophone  interval  8-9  (labeled  8).  If  we  approximate  this 
curve  as  a  straight  line,  it  increases  from  to  3500  m/s  in  about  44  hours.  In  Fig.  6  the  slope  of 
the  curve  between  16.9  and  22.0  cm  (the  depths  of  hydrophones  8  and  9)  gives  a  freezing  rate  of 
0.15-0.16  cm,/h.  The  speed  starts  deviating  from  as  ice  begins  to  form  in  the  interval  between 
hydrophones.  For  an  initial  estimate,  we  assume  that  the  skeletal  layer  has  completely  traversed  the 
interval  when  the  speed  reaches  3500  m/s.  With  this  assumption  and  the  spacing  between  hydro¬ 
phones  from  Fig.  2,  curve  8  gives  a  skeletal  layer  thickness  of  1.5  to  1.9  cm.  Similarly,  cur\e  5 
gives  a  thickness  of  2.7  to  3. 1  cm.  These  are  clearly  underestimates  of  the  layer  thickness  since  the 
lower  hydrophone  in  an  interval  would  still  be  in  the  skeletal  zone  when  the  average  speed  in  the 
interval  reaches  3500  m/s.  Furthermore,  the  underestimate  is  more  severe  for  curve  8  since  the  in- 
ters-al  between  hydrophones  8  and  9  is  about  twice  that  between  h>drophones  5  and  6.  The  next 
subsection  uses  the  entire  datti  set  to  obtain  both  ;•  better  estimate  of  the  skeletal  layer  thickness  and 
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FIG  9.  The  change  in  sound  speed  as  measured  for  each  hydrophone  p; 
ering  one  hydrophone  after  another.  The  number  of  the  deeper  of  the  pa 
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an  estimate  of  the  velocity  profile  within  the  skeletal  layer. 

4.  Sound  speed  in  the  skeletal  layer 

Every  set  of  sound  speed  measurements  included  a  measurement  across  the  lowest  few  centi¬ 
meters  of  ice.  We  assumed  the  skeletal  layer  formed  consistently  at  all  depths,  and  combined  all 
data  to  obtain  a  sound  speed  profile  for  the  layer. 

The  sound  speed  C/(5)  in  ice  as  a  function  of  distance  (s)  away  from  the  ice/water  interface  was 
determined  by  first  considering  hydrophone  intervals  that  included  the  interface  (Fig.  10a).  To  de¬ 
termine  the  extent  of  the  ice  below  a  hydrophone,  the  depth  of  the  interface  at  the  measurement 
time  was  taken  from  Fig.  6.  An  average  value  for  the  acoustic  speed  Cjis)  in  this  portion  of  ice  was 
then  calculated  using  and  the  average  speed  Q  in  the  interval  between  the  adjacent  hydro¬ 
phones.  To  get  larger  values  of  s,  the  interval  was  extended  to  other  hydrophones  above  (Fig.  10b); 
the  average  speed  was  determined  by  adding  travel  times  for  intervals  up  to  15  cm.  These  are  plot¬ 
ted  as  points  in  Fig.  10c.  These  were  fitted  with  a  fourth  order  polynomial  (dotted  line  in  Fig.  10c) 
giving,  for  s  in  centimeters, 

Cjis)  =  1490  +  480.93.? -  55.67(h- .3. 254S.v-().()7 1579. 

which  is  a  good  fit  for  an  s  of  0  to  15  cm.  The  total  travel  time  through  the  ice  can  be  expressed  as 
either  an  integral  or  in  terms  of  the  average  speed  given  in  Eq.  (1).  Taking  the  derivative  of  this 
travel  time  equality  with  respect  to  .v  gives 

d  I  [ !  d  i  ■'f  'I 

~  j  (2) 
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or 

—  2 
C,  (5) 

^  {C,(s) -s^jC,(s))' 

which  is  the  result  plotted  as  a  solid  line  in  Fig.  10c.  It  is  shown  only  to  .v=7  cm;  beyond  this  point 
it  becomes  erratic  because  we  have  only  sparse  data  for  the  large  s  values  in  the  fourth  order  fit  of 
Eq.(l). 

To  obtain  sound  speeds  farther  into  the  ice,  average  speeds  for  some  hydrophone  pairs  were 
determined  in  intervals  where  the  speed  had  reached  a  plateau.  The  values  obtained  are  plotted  as 
open  circles  in  Fig.  10c  and  extend  the  speed  profile  beyond  the  end  of  the  solid  line.  The  result  is 
a  high  resolution,  in-situ,  vertical  profile  for  the  speed  of  the  longitudinal  wave  in  growing  arctic 
sea  ice. 

The  “acoustic”  skeletal  layer  thickness  apparent  from  the  solid  line  is  about  3  cm.  As  expected, 
this  is  larger  than  our  initial  estimates  discussed  earlier.  We  use  the  term  “acoustic”  skeletal  layer 
because  the  thickness  measured  here  may  differ  from  that  obtained  by  examining  the  ice  structure 
directly.  However,  the  value  found  is  commensurate  with  those  reponed  in  Ref.  13  (3  cm)  and  Ref. 
14  (1.8  to  5.75  cm), 

5.  Absorption  of  sound 

The  hydrophone  data  were  also  used  to  calculate  the  absorption  of  sound  in  the  ice  at  different 
temperatures.  Recording  each  ping  on  two  adjacent  hydrophones  simultaneously  allowed  a  com¬ 
parison  of  amplitudes  with  no  effects  from  ping-to-ping  variability  in  the  transmitter.  The  relative 
sensitivities  of  the  hydrophones  in  water  were  obtained  by  recording  signals  before  any  ice  had 
formed.  We  calculated  the  absorption  of  sound  in  the  ice  assuming  that  the  .sensitivity  of  a  hydro¬ 
phone  did  not  decrease  when  ice  formed  around  it.  If  both  hydrophones  were  in  ice,  both  sensiiis  - 
ities  would  change  about  the  same,  and  our  calculated  absorption  would  be  accurate.  Howes  er.  if 
only  the  upper  hydrophone  was  in  ice  and  its  sensitivitx  did  indeed  decrease,  our  calculated  absorp¬ 
tion  for  the  bracketed  ice,  which  generally  includes  ilie  skeletal  laser,  would  be  high 


FIG.  10.  The  average  sound  speed  for  a  distance  s  into  the  ice  is  calculated  using  the  measured  trav¬ 
el  times  betv^een  a  hydrophone  just  below  the  ice  surface  and  one  just  above  as  in  (a).  To  obtain 
larger  5  values,  the  second  or  third  hydrophone  above  the  ice  surface  w  as  substituted  for  the  upper 
one.  as  in  (b;.  This  gives  two  or  three  values  of  average  sound  speed  \  ersus  5  for  every  set  of  mea¬ 
surements  during  the  ice  grow  th  experiment,  these  are  plotted  as  dots  in  (c).  From  a  least  squares 
fit  of  these  data.  Eq.  (1 1.  the  profile  in  the  ice  has  been  calculated.  Eq  and  shown  as  a  solid 
line.  The  line  is  extended  to  fit  some  sound  speed  measurement.^  (circlesi  obtained  from  hydro¬ 
phone  pairs  well  into  the  ice 
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For  adjacent  hydrophones  a  distance  (s)  centimeters  apart  with  amplitude  readings  of  A'  and 
B'  in  water,  and  A  and  B  after  partial  or  complete  embedment,  the  absorption  loss  (in  decibels)  is 
calculated  as 


AB' 

L,  =  201oi.{--J. 


(4) 


This  total  absorption  loss  can  be  expressed  in  terms  of  an  interface  loss  p  and  an  absorption  coef¬ 
ficient  per  centimeter  a  as 


S 

=  p  +  ja(s)ds 
0 


(5) 


If  we  use  the  temperature  dependence  of  the  McCammon-McDaniel  equation^  for  the  absorption 
coefficient,  i.e.. 


(S  ^ 

a  =  i/i-p 


(6) 


(where  T  is  temperature  in  degrees  Celsius,  k  isu  constant,  and  /is  frequency  in  kilohertz)  and  as¬ 
sume  a  linear  temperature  profile  in  the  ice,  Eq.  (5)  can  be  written  as  an  integral  over  temperature 
to  obtain 


Z.,  =  P^kfqs, 


(7) 


where  q  is  the  integrated  temperature  effect  for  thickness  .v.  Witlt  temperature  T2  at  the  upper  hy¬ 
drophone  and  T /  at  the  lower  hydrophone. 


(S) 
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In  Eq.  (6)  intergrain  friction  is  assumed  to  be  the  cause  of  the  acoustic  absorption.  Again,  this 
is  most  probable  away  from  the  skeletal  layer.  Within  the  layer,  relative  fluid/solid  motion  (as  ac¬ 
counted  for  in  the  Biot  model  of  a  porous  medium)  may  dominate'^  over  frame  loss,  such  as  inter- 
grain  friction,  giving  large  attenuations  with  a  different  frequency  dependence  than  in  Eq.  (6). 
Equation  (7)  shows  that  a  plot  of  measured  loss  vs  fqs  should  display  a  linear  variation  with  slope 
k.  If  temperatures  along  the  path  are  measured,  q  can  be  determined  more  accurately  by  numerical¬ 
ly  integrating  over  each  hydrophone  interval. 

Pairs  of  amplitude  readings  were  taken  for  each  interval  before  any  ice  appeared  in  the  interval. 

Their  ratios  were  fairly  consistent,  showing  that  the  equipment  was  capable  of  repeatable  loss  mea¬ 
surements.  Averages  of  these  ratios  over  several  measurements  provide  a  calibration  of  the  receiv¬ 
ing  circuits. 

We  now  examine  the  absorption  loss  as  the  freezing  progressed.  Data  for  hydrophones  10-12 
are  unusab.e,  as  mentioned  earlier.  Proceeding  downward,  losses  were  measured  for  all  hydro¬ 
phone  intervals  in  which  ice  was  present.  Figure  1 1  shows  how  the  loss  for  travel  from  hydrophone 
7  to  8  and  8  to  9  varied  throughout  the  severaj^runs^^seg=tgp  of  Fi^  JrfnrTimr-ofatttA  A  plot  of  the 
sound  speed  is  shown  for  comparison.  The  high  absorption  appears  to  occur  where  the  sound  speed 
is  increasing,  i.e.,  within  the  skeletal  layer.  (As  a  point  of  reference,  measurements'^  we  have  made 
by  other  techniques  indicated  a  total  absorption  in  the  skeletal  layer  of  about  2.3  dB.) 

Absorption  in  the  skeletal  layer  can  be  examined  using  the  entire  data  set.  In  every  set  of  mea¬ 
surements  the  lower  surface  of  the  ice  lies  between  two  hydrophones.  Thus,  in  this  hydrophone  in¬ 
terval  there  is  both  ice  and  water.  The  temperature  at  the  hydrophone  above  the  ice  interface  is 
interpolated,  assuming  a  constant  gradient  between  thermistors.  This  and  the  water  temperature  are 
used  to  compute  q  for  the  layer.  The  transmission  loss  between  the  hydrophones  is  assumed  to  oc¬ 
cur  in  the  ice,  of  thickness  s.  When  the  loss  is  computed  for  many  measurements,  and  thus  many 
values  of  s,  we  have  a  relationship  showing  the  loss  in  various  thicknesses  of  the  skeletal  layer. 

Some  additional  data  points  were  obtained  by  taking  an  intercal  involving  tuo  hydrophones,  sum¬ 
ming  the  two  losses,  and  thus  obtaining  some  larger  values  for  .v.  All  data  are  plotted  in  Fig  12 
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FIG.  1 1 .  Signal  loss  and  sound  speed  in  the  interv-al  between  hydrophones  pairs;  (a)  7-8  and  (b)  8- 
9.  The  arrow  s  indicate  when  each  hydrophone  was  frozen  into  the  ice.  When  a  hydrophone  pair  is 
sampling  the  skeletal  layer,  the  loss  is  large — perhaps  excessive — since  it  includes  effects  due  to 
sensitivity  change  As,  the  ice  advances,  the  hydrophones  are  left  in  solid  ice  which  has  a  low  ab¬ 
sorption  The  sound  speed  gradualls  increases  as  more  and  more  of  the  interv  al  becomes  ice. 
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A  best-fit  second  order  curve  for  the  data  shows,  by  its  intercept,  a  boundary  transmission  loss 
of  4.5  dB.  Since  the  speed  is  faster  in  the  ice,  implying  that  the  pressure  transmission  coefficient 
should  be  greater  than  1  (assuming  no  density  change),  this  loss  must  be  primarily  an  artifact  of 
reduced  hydrophone  sensitivity.  The  slope  of  the  curve  near  qs=0  gives  a  peak  loss  rate  of 
fc9’dB/cm,  an  upper  limit  to  the  skeletal  layer  absorption  coefficient,  r?  shows  a  total  ab- 


sorption  in  the  0-3  cm  skeletal  layer  of  4  dB,  is  a  little  higher  than  the  2.3  dB  quoted  above  from 

A 

previous  work.  A  fourth  degree  curve  has  been  added  (dotted  line)  to  provide  a  better  fit  in  the  skel- 

etal  layer.  This  curve  indicates  a  peak  loss  of  »b^'3=^dB/cm  in  the  skeletal  layer.  A  statistical  F- 

A-  >1  yi 

test  of  the  data  in  Fig.  12  indicates  that  the  decrease  in  X  using  the  second-order  fit  (as  compared 
to  a  linear  fit)  is  significant  at  the  94th  percent  level.  For  the  fourth-order  fit,  the  decrease  is  signif¬ 
icant  at  the  9()th  percent  level. 

The  slope  of  the  curve  at  any  point  also  detemiines  the  value  of  the  coefficient  k  mentioned 
above,  which  goes  from  1 .0  (2.5  for  the  fourth  degree  curve)  near  the  ice/water  interface  to  0.3.  We 
reiterate  that  k  would  not  vary  if  the  form  of  the  absorption  throughout  the  ice,  including  the  skel¬ 
etal  layer,  was  as  given  in  Eq.  (6).  Taken  together,  our  previous  work  and  the  data  in  Figs.  1 1  and 
12  imply  a  failure  of  this  temperature  dependence  in  the  skeletal  layer. 
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6.  Reflections  from  the  underside  of  ice 

Interleaved  between  the  92  kHz  transmission  measurements  from  the  phnter  transducer  to  the 
hydrophones  were  measurements  of  the  return  to  the  platter  from  the  underice  surface  at  several 
frequencies.  A  comparison  with  initial  water  surface  retu.ns  gave  an  amplitude  reflection  coeffi¬ 
cient. 

Coefficients  were  calculated  for  each  set  of  measurements.  Because  of  rather  high  variability, 
an  average  was  taken  for  the  several  runs  of  each  day.  These  are  plott.  J  in  Fig.  13.  The  220  kHz 
values  were  too  variable  and  have  been  omitted  (the  beam  is  very  narrow  ai  220  kHz,  and  the  spot 
size  may  have  been  too  small  to  provide  adequate  averaging). 

Figure  14a  gives  averages  of  the  reflection  coefficicnis  measincd  during  freezing,  along  uiih 
results  from  earlier  measurements  for  the  undersurface  i>r  the  cano]y\  The  coefficients,  and  tlie  de- 
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FIG.  1 2.  Signal  loss  measured  between  hydrophones  and  calculated  for  various  distances  s  into  the 
ice.  The  loss  is  plotted  versus  the  integral  of fqs.  where /  is  the  frequency  and  q  includes  the  tem¬ 
perature  dependence  of  the  absorption.  The  slope  in  this  plot  is  the  absorption  constant  k  in  Eq.  (6). 
A  second  degree  best  fit  curve  shows  that  k  is  high  in  the  skeletal  layer  and  then  drops  off  fanher 
into  the  ice.  A  founh  degree  curve  (.dotted)  gives  a  better  fit  for  the  ?-cm  skeletal  layer. 
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FIG.  13.  Daily  averages  of  the  amplitude  reflection  coefficients  at  several  frequencies  obtained 
during  ice  growth 
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crease  with  frequency,  agree  well.  Also  shown  is  the  theoretical  prediction  of  the  reflection  coef¬ 
ficient  (labeled  if)  using  the  sound  speed  profile 

A 

-f  ‘  r. 

I  C4-H  (C3-C4)  { 1 -lanh  ('6l(5-jrQV] }  +  ^(^2-^4)  { 1  +  lanh  1 


with 

^0  =  0.0  m,  =  0.065  m,  C)  =  1436  m/s,  Co  =  4000  m/s,  C3  =  1490  m/s. 


and  a  density  of  0.959  g/cm^.  Equation  (9)  is  plotted  in  Fig.  14b,  with  the  curve  from  Fig.  10c 
shown  as  a  dotted  line  for  comparison.  Here,  Cj  is  sound  speed  in  water,  C2  is  bulk  sound  speed  in 
the  interior  of  the  ice,  cj  is  sound  speed  in  the  ice  at  the  ice/water  interface,  and  C4  is  sound  speed 
at  one  point  between  .Vg  and  The  theoretical  calculation^^  is  discussed  in  Section  IV  of  Ref.  8 
where  it  was  used  to  predict  the  sound  speed  profile  in  the  ice  using  the  reflection  coefficient  data 
for  1988  shown  in  Fig.  14.  The  skeletal  layer  was  predicted  to  be  2  to  3  cm  thick  (Fig.  16a  of  Ref. 
8).  The  comparison  shown  here  between  reflection  data  and  theoretical  calculation  using  a  known 
sound  speed  profile  demonstrates  the  utility  of  the  type  of  acoustic  inversion  performed  in  Ref.  8. 

B.  Submerged  Ice  Block 
1,  Ice  block  appearance 

The  second  phase  of  the  experiment  was  initiated  by  coring  out  a  58-cm  diameter  cylindrical 
block,  centered  on  the  sensor  arrays,  using  an  APL  thermal  drill, The  block  was  then  held  down 
with  a  rigid  structure.  An  underwater  videotape  of  the  block,  taken  after  the  experiment  was  com¬ 
pleted,  clearly  showed  that  all  surfaces  of  the  block  were  fuzzy  due  to  ice  growth.  Efforts  to  repro¬ 
duce  the  view  for  this  article  resulted  in  Fig.  15,  in  which  everything  looks  fuzzy.  Above  the  block 
is  the  inpod  structure;  the  ends  of  the  sensor  tirrays  are  visible  at  the  bottom  of  the  block. 


l.S 


FIG.  14.  In  (a),  the  average  amplitude  reflection  coefficients  for  all  measurements  during  ice 
growth  (filled  circles)  are  compared  with  ice  canopy  reflectivities  measured  in  1988.  In  (b),  the 
dotted  line  is  the  sound  speed  profile  derived  from  measurements  (the  heavy  curve  of  Fig.  10c). 
The  solid  line  in  (b)  was  calculated  using  Eq.  (9)  with  parameters  selected  to  obtain  a  good  repre¬ 
sentation  of  the  data;  it  was  used  to  predict  the  curve  in  (a). 


FIG.  15.  An  undervk-ater  viev^  of  the  submerged  ice  block. 
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2.  Temperature  monitoring  and  calculation  of  thermal  diffusivity 

The  readings  of  temperature  in  the  ice  block  were  continued  after  it  was  submerged.  The  pro¬ 
files  obtained  are  plotted  in  Fig.  16.  On  the  last  set  of  readings  the  entire  block  was  near  the  tem¬ 
perature  of  the  water.  We  used  profiles  taken  0.4,  8,  and  16  hours  after  submergence  (run  numbers 
1,  17,  and  19)  to  calculate  an  “effective”  thermal  diffusivity  (see  appendix),  assuming  no  spatial  or 
temporal  variance  in  the  block’s  thermal  diffusivity.  The  value  obtained  is  between  0.0013  and 
0.0018  cm^/s,  about  one-fifth  the  value  calculated  earlier. 

A  two-dimensional  array  of  thermistors  would  allow  a  more  complete  calculation  of  the  diffu¬ 
sivity,  taking  into  account  the  temporal  variability  of  the  thermal  diffusivity  within  the  block.  How¬ 
ever,  one  can  qualitatively  explain  the  low  value  of  effective  diffusivity  in  the  submerged  block  by 
examining  the  equation  given  for  the  dependence  of  diffusivity  on  salinity,  temperature,  and  den¬ 
sity  [Eq.  (A  16)  in  the  appendix]  and  using  an  electrical  analogy. For  the  present  temperatures 
(-10°C  to  -1.8°C)  the  term  in  brackets  in  the  denominator  (the  specific  heat)  varies  most  dramati¬ 
cally.  It  is  about  ten  times  larger  in  the  outer  portion  of  the  ice,  implying  a  much  lower  diffusivity 
there.  Using  the  electrical  analogy  of  Ref.  19  and  treating  the  block  in  layers,  one  finds  that  the 
capacitance  of  the  outer  portion  of  the  block  is  in  parallel  with  the  inner  portion,  and  the  time  rate 
of  change  of  the  temperature  is  dominated  by  the  high  capacitance  region.  Physically,  the  outer 
portion  of  the  block  can  absorb  a  large  amount  of  heat  via  pha.se  changes  (latent  heat  is  part  of  the 
specific  heat  given  here)  with  a  .small  change  in  temperature.  This  tends  to  shield  the  inner  ponion. 
Thus  the  temperature  change  there  is  small.  More  fomially,  the  high  capacitance  in  the  outer  por¬ 
tion  affects  the  inner  portion  temperature  change  due  to  the  “boundary”  conditions  of  continuity  of 
temperature  and  heat  conduction  between  the  inner  and  outer  portions  of  the  ice. 

3.  Sound  speed  monitoring 

The  submergence  phase  of  the  experiment  was  divided  into  seven  time  periods.  For  each,  a 
sound  speed  was  computed  in  each  hydrophone  interval  by  averaging  the  data  obtained  in  that  pe¬ 
riod.  The  results  are  plotted  in  Fig.  17.  At  most  depths  there  was  a  rapid  change  in  sound  speed 
followed  by  a  slow  drift  as  the  block  was  warmed  by  the  water.  The  souiui  speeds  for  the  deepest 
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FIG  16.  Temperature  profiles  in  the  submerged  ice  block  during  the  four  days  it  took  for  the  block 
to  w.arm  to  the  temperature  of  the  lAater.  The  four  profiles  taken  before  coring  are  the  same  four 
shou.n  as  dotted  in  Fig.  8.  As  in  Fig  8.  several  profiles  are  labeled  b>  run  number,  and  the  graph 
at  the  top  of  the  figure  associates  run  number  with  local  time 
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FJG  ]7.  Sound  speed  between  hsdrophone  pairs.  obt;uned  b\  averaging  over  successive  time  in¬ 
tervals.  The  first  measurement  took  place  2  hours  and  20  min  after  the  block  was  submerged. 
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hydrophone  interval  increase  with  time  because  additional  ice  fomied  in  the  interval  after  submer¬ 
gence  owing  to  heat  conduction  into  the  ice  above.  At  the  upper  levels  the  sound  speed  decreases 
as  the  ice  warms.  For  hydrophone  pairs  3-4, 4-5,  and  5-6  the  speed  goes  through  maxima  after  sub¬ 
mergence. 

An  alternate  view  of  the  hydrophone  data  is  shown  in  Fig.  18.  Here  sound  speed  profiles  in  the 
ice  before  submergence  and  for  several  time  intervals  after  submergence  are  derived  from  the  data. 
Note  that  the  sound  speed  maxima  for  the  first  day  in  hydrophone  pairs  3-4,  4-5,  and  5-6  are  also 
associated  with  a  maximum  in  the  sound  speed  profile.  Indeed,  one  sees  that  from  the  bottom  of 
the  block  until  at  least  15  cm  above  this  interface  there  is  an  increase  of  speed  with  time  during  part 
of  the  experiment,  the  greatest  increase  occurring  about  1 1  cm  from  the  bottom.  The  data  again 
indicate  that  some  shon  time  scale  effect  causes  the  profile  to  have  a  maximum  which  then  relaxes 
via  a  much  longer  time  scale  mechanism. 

The  long  time  scale  relaxation  of  the  profile  can  be  attributed  to  warming  of  the  block.  The 
short  time  scale  effect  was  unexpected.  At  present  we  believe  it  can  be  explained  by  disruption  of 
the  hydrostatic  equilibrium  when  the  block  is  pushed  down  into  the  water.  To  support  this  hypoth¬ 
esis,  we  first  recall  the  analysis  of  Hide  and  Manin,^"^  who  have  calculated  the  pressure  within  a 
brine  drainage  channel  and  the  hydrostatic  pressure  at  a  height  z  above  the  ice/water  interface  to 
explain  the  expulsion  of  brine  in  young  growing  sea  ice.  Using  their  Eq.  7a  and  b,  the  volume  flux 
iQ)  out  of  a  model  brine  drainage  channel  of  constant  radius  (r)  is 


Q 


nr 


(10) 


where  is  the  mean  density  of  the  brine  in  the  channel,  v^^  is  the  mean  viscosity  of  the  brine, 
is  the  level  of  the  brine  in  the  channel,  is  the  brine  hydrostatic  pressure  at  the  bottom  of  the 
channel,  and  p  is  the  .seawater  hydrostatic  pressure  at  the  bottom  of  the  channel.  During  the  freez¬ 
ing  process,  gravity  drainage  of  the  heavier  brine  to  obtain  pressure  equilibrium  at  the  lower  sur¬ 
face  of  the  ice  results  in  voids  in  the  upper  pan  of  the  drainage  channeK  In  submerging  the  block 
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FIG  1 8.  Sound  speed  measurements  for  the  submerged  phase.  The  profile  labeled  0  is  the  last  one 
taken  before  coring  of  the  block;  profile  1  is  the  average  of  four  measurements  5  h  after  submer¬ 
gence;  2  is  after  12  h  of  warming;  3  is  after  1  day’s  warming;  and  4  was  taken  after4daysof  warm¬ 
ing. 
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such  that  the  top  is  at  1  m  depth,  we  have  increased  p  in  Eq.  (10)  by  about  10^  N/m^  (1.5  Ib/in^). 
With  the  assumption  that  the  upper  surface  of  the  ice  is  impermeable  compared  with  the  lower  part 
of  the  ice,  this  results  in  a  flux  of  low  salinity  sea  water  into  the  brine  channel,  lowering  the  salinity 
of  the  brine.  Because  this  lower  salinity  brine  is  now  out  of  thermal  equilibrium,  freezing  occurs 
and  brings  the  brine  salinity  up  to  its  equilibrium  value. ^  The  hypothesis  is  that  this  e  Tect  produces 
the  increase  in  speed  seen  in  the  bottom  part  of  the  ice,  in  a  manner  to  be  discussed  below.  If  this 
is  the  correct  explanation,  the  region  near  the  lower  water/ice  interface  (15  cm  thick)  that  experi¬ 
ences  an  increase  in  sound  speed  may  be  the  region  where  brine  drainage  channels  occupy  a  sig¬ 
nificant  volume  of  the  bulk  ice.  That  distance  is  about  one-third  of  the  ice  thickness  in  this  case. 

Porosity  is  one  of  the  major  parameters  determining  the  speed  of  sound  in  ice.®’*  ^  One  mech¬ 
anism  by  which  a  maximum  could  be  developed  in  the  vertical  profile  is,  therefore,  if  the  salinity 
reduction  by  seawater  influx  was  sufficient  to  reduce  the  porosity  by  the  necessary  amount.  At  the 
end  of  the  submergence  phase  we  took  an  ice  core  and  measured  the  salinity  profile.  Figure  19 
shows  the  profiles  obtained  then  and  at  the  end  of  the  ice  growth  phase.  These  profiles  and  the  ap¬ 
propriate  temperature  profiles  were  used  to  obtain  the  porosity  profiles  shown  in  the  figure.  Also 
displayed  is  sound  speed  profile  No.  4  from  Fig.  18.  The  lower  salinity  near  the  peak  in  sound 
speed  at  30  cm  depth  qualitatively  confirms  the  physical  picture  given  above.  However,  the  effect 
on  porosity  is  not  sufficient  to  predict  the  speed  enhancement  using  present  model  parameters.  To 
demonstrate  this,  in  Fig.  20  we  compare  the  experimental  results  with  the  results  using  the  poros- 
ities  of  Fig.  19  in  two  different  models  of  porosity  vs  sound  speed.  ’  Figure  20a  shows  model  vs 
experiment  at  the  end  of  the  ice  growth  phase;  Fig.  20b  shows  model  vs  experiment  at  the  end  of 
the  submergence  phase.  Only  a  slight  increase  is  seen  in  the  model  results  for  the  area  of  reduced 
salinity;  a  large  increase  is  obserx'ed  in  the  experimental  results. 

The  porous  mode!  of  ice  in  Ref.  1 1  suggests  another  mechanism  for  the  increase  in  speed.  One 
of  the  parameters  in  that  model,  in  addition  to  porosity,  is  the  bulk  modulus  of  the  solid  ice  frame. 
If  the  freezing  process  causes  a  more  rigid  frame  structure,  the  acoustic  speed  increases  Since  the 
freezing  that  would  occur  via  the  inllux  of  salt  water  upon  submergence  i^  cenainls  not  tlie  ts  pical 


Ih 


Cm/s)  Cppi  )  C®/o) 


FIG.  19.  Measured  salinities  of  the  ice  cores  before  submergence  (SI)  and  after  warming  (S2)  were 
used  to  calculate  porosit\ .  The  sound  speed  profile  from  Fig.  18  (No.  4)  is  plotted  to  show  how  the 
peak  at  32  cm  corresponds  to  a  low  salinity. 
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FIG.  20.  Measured  sound  speed  profiles  (a)  before  and  (b)  after  4-day  submersion  are  compared 
with  those  calculated  from  measured  temperature  and  salinity  profiles  in  ice  cores  using  the  models 
in  Ref.  8  (dashed  line)  and  Ref.  1 1  (doned  line).  Note  that  the  reduced  salinity  seen  in  Fig.  19  is 
not  sufficient  to  produce  the  high  sound  speed  seen  near  the  bonom  of  the  ice  after  submersion  for 
4  days. 
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freezing  process  in  ice  growth,  it  is  possible  that  the  frame  rigidity  is  different.  There  is  some  pre¬ 
vious  evidence  of  variations  in  speed  not  related  to  porosity.  In  Ref.  20  a  different  Young’s  moduli 
versus  porosity  curve  was  found  for  winter  ice  than  for  summer  ice.  The  authors  commented  that 
the  Young’s  modulus  might  also  depend  on  temperature.  Within  the  porous  ice  model  this  is  equiv¬ 
alent  to  saying  the  frame  rigidity  depends  on  temperature,  a  contention  that  seems  reasonable. 

4.  Absorption  changes 

As  in  the  ice  growth  phase  of  the  experiment,  sound  pulses  transmitted  into  the  block  from  be¬ 
low  were  also  used  to  determine  absorption.  For  each  ping  we  compared  the  amplitudes  at  two  ad¬ 
jacent  hydrophones.  The  ratio  was  compared  with  the  in-water  ratio  obtained  for  these  two 
hydrophones  before  the  ice  began  to  freeze. 

Temperatures  were  measured  every  few  hours.  There  were  twice  as  many  thermistors  as  hy¬ 
drophones,  so  the  temperature  gradient  between  hydrophones  was  well  known.  We  calculated  the 
absorption  loss  from  the  amplitude  ratio  as  in  Eq.  (4).  For  each  run,  a  plot  of  cumulative  absorption 
loss  vs  cumulative  qs  for  hydrophones  1-9  gave  eight  points.  A  best  fit  line  gave  a  slope  equal  to 
k.  The  average  k  was  (0.3  ±  0.1).  This  is  60%  higher  than  our  1988  results,'^’  but  in  good  agreement 
with  that  shown  in  Fig.  12.  The  variation  of  k  during  the  5  days  of  submergence  is  shown  in  Fig. 
21.  Individual  k's  are  erratic.  There  is  no  appreciable  trend  in  k  with  time,  which  suggests  that  the 
temperature  dependence  given  by  McCammon  and  McDaniel  [Eq.  (6)|  is  appropriate.  This  test  is 
a  moderately  good  one  because  the  dependence  gives  a  doubling  of  the  absorption  for  the  5- 
day  temperature  change  from  an  average  of  -6”C  to  -2°C. 

5.  Underice  reflectivity  changes 

After  the  block  was  cored  out  and  depressed,  readings  of  the  return  from  the  lower  surface  of 
the  block  were  obtained  to  monitor  the  effect  of  block  warming  on  the  amplitude  of  the  reflection 
coefficient.  The  return  from  the  block  w'as  converted  to  target  strength  and  tiien  compared  with  the 
theoretical  target  strength  of  a  5X-cm  diameter  rigid  plate.  The  difference  is  20  log  Ra.  where  Af; 
is  the  amplitude  refection  coefficient.  It  is  computed  and  plotted  in  Fig.  22,  wiiich  show  s  a  decrease 
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FIG  2 1  Measurements  of  absorption  coefficient  k  dunng  the  4  days  of  block  submersion.  A  least 
squares  line  shov.^,  only  a  small  slope,  supponing  the  temperature  dependence  of  Eq.  (6). 


FIG.  22.  Measured  amplitude  reflection  coefficients  during  the  wanning  phase.  Although  the  scat¬ 
ter  of  individual  measurements  was  large,  averages  for  each  day  show  a  decreasing  trend  after  the 
initial  increase  when  the  block  is  submerged. 
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in  Ra  as  the  block  warms.  Figure  23a  gives  the  reflection  coefficient  as  a  function  of  frequency  at 
the  end  of  the  experiment,  compared  again  with  the  1988  results  and  experimental  results  in  Fig. 
14a.  Also  shown  is  a  theoretical  prediction  calculated  using  Eq.  (9)  but  with 

^0  =  0.0  m,  5,  =  0.16  m,  Ci  =  1436  m/s,  C2  =  3300  m/s,  03=  1500  m/s, 

C4  =  4400  m/s,  5o  (C4  -  C3)  =  50. 000  (1/s),  V2  5,  (C2  -  C4)  =  -15000  (1/s). 

Equation  (9)  with  the  parameters  above  is  plotted  in  Fig.  23b  and  compared  with  the  experimental 
profile  shown  in  Fig.  20b.  Again,  data  and  theory  agree  well;  both  show  a  decrease  due  to  the  re¬ 
laxation  of  the  speed  gradient  at  the  bottom  of  the  ice. 

III.  SUMMARY 

We  have  presented  an  overview  of  some  of  the  in-sitii  ice  properties  of  growing  and  newly  sub¬ 
merged  sea  ice.  Together,  the  ice  growth  and  subsequent  block  submergence  approximate  the  life 
cycle  of  an  individual  block  in  a  newly  fomied  ice  keel. 

One  of  the  main  motivations  of  this  work  was  to  examine  the  acoustic  properties  of  highly  po¬ 
rous  regions  of  ice.  In  this  regard,  the  high  resolution  sound  speed  profiles  shown  in  Figs.  10  and 
18  are  of  primary  importance.  Such  in-situ  results  are  needed  for  predicting  acoustic  scattering 
from  the  arctic  canopy-a  specific  example  is  the  prediction  of  nomial  incidence  reflectivity  as 
shown  in  Figs.  14  and  23.  The  results  for  the  growing  ice  (Fig.  10)  show  an  acoustic  skeletal  layer 
(defined  as  having  an  acoustic  speed  of  <3500  m/s)  with  a  thickness  of  about  3  cm.  The  concurrent 
measurement  of  scattering  from  the  growing  ice,  to  determine  reflection  coefficients,  demonstrated 
the  viability  of  sensing  remotely  from  below  to  determine  the  longitudinal  wave  acoustic  speed  of 
the  sea  ice.  The  frequency  dependence  of  the  reflection  coefficient  measured  during  the  ice  growth 
phase  agrees  with  previous  work,  and  the  sound  speed  profiles  inferred  from  that  work  are  very 
close  to  the  ones  shown  here  (Fig.  10). 

Absorption  was  also  measured  in  situ.  The  data  show  peak  values  in  the  skeletal  layer  between 
2  and  5  dB/cm  at  92  kHz.  Using  the  temperature  dependence  of  the  McCammon-McDaniel  equa¬ 
tion  [Eq.  (6)1,  the  constant  k  decreases  from  a  high,  poorly  defined  value  (1-3)  at  the  interface  to, 
for  example,  a  value  of  0.3  at  20  cm  into  the  ice.  In  the  submerged  ice  block,  a  k  of  0.3  w  as  deter- 
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FIG.  23.  In  (a),  the  average  amplitude  reflection  coefficients  for  all  measurements  during  the  sub¬ 
merged  phase  (filled  circles)  are  compared  with  ice  canopy  reflectivities  measured  in  1988.  The 
dotted  line  represents  the  data  from  Fig.  14  for  the  ice-grov*  th  phase.  In  (b),  the  sound  speed  data 
for  the  end  of  the  submerged  phase  is  a  re-plot  of  sound  speeds  shown  by  the  #4  curve -in  Fig.  18. 
The  line  is  a  fit  of  Eq.  (9)  to  the  data,  and  this  profile  was  used  to  predict  the  solid  line  in  (a). 
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mined  as  the  block  warmed  to  the  water  temperature,  thus  providing  a  check  on  the  temperature 
dependence  seen  in  Eq.  (6).  The  scatter  of  the  data  for  k  during  this  period  (Fig.  21)  is  large  but  no 
noticeable  trend  is  seen.  This  indicates  that  the  temperature  dependence,  which  is  a  result  of  as¬ 
suming  that  intergra  .i  interaction  controls  the  absorption,  is  useful  in  this  case.  This  implies  that, 
though  the  porosity  increases,  a  pore  structure  remains  which  prevents  other  effects  (such  as  fluid 
viscosity  in  the  Biot  model)  from  altering  the  temperature  dependence  of  the  absorption. 

The  measured  temperature  profiles  in  the  ice  were  used  to  calculate  the  thermal  diffusivity  in 
the  upper  portion  of  growing  ice  (0.0080  cm^/s)  and  in  the  submerged  ice  block  (about 
0.0015  cm  /s).  The  difference  between  these  values  was  related  to  the  change  of  specific  heat  of 
the  ice  as  a  function  of  temperature.  Measured  salinity  and  temperature  profiles  allowed  calcula¬ 
tion  of  porosities  in  the  growing  and  submerged  ice.  These  porosities  were  then  used  in  two  mod¬ 
els.  Figure  20a  shows  comparisons  between  model  and  experiment  during  ice  growth.  The 
subsequent  comparison  to  the  experimental  profile  in  the  submerged  block  (Fig.  20b)  suggests  that 
porosity  changes  alone  are  not  responsible  for  the  high  speed  regime  seen  there. 

Sound  speed  profiles  taken  during  the  submergence  ponion  of  the  experiment  (Fig.  18)  show- 
effects  from  long  (up  to  80  hours)  and  short  (a  few  hours)  time  scale  processes.  The  long  time  scale 
effect  corresponds  to  the  time  scale  of  block  warming.  The  short  time  scale  process  is  hypothesized 
to  be  the  result  of  disruption  of  the  hydrostatic  equilibrium  by  block  submergence.  We  believe  the 
submergence  forces  seawater  into  the  block,  thereby  reducing  the  salinity  in  brine  drainage  chan¬ 
nels  within  the  bottom  third  of  the  ice,  which  reduces  the  porosity  through  freezing  and  also  chang¬ 
es  the  ice  frame  rigidity.  A  reduced  salinity  in  the  region  of  increased  sound  speed  supports  this 
hypothesis.  Additional  information  regarding  this  work  is  given  in  Ref  23. 
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APPENDIX  -  Calculations  of  Thermal  Diffusivity 
During  Ice  Growth 

We  use  the  warming  trend  seen  in  Fig.  4  (29  March)  and  the  resulting  change  in  temperature 
profile  shown  in  Fig.^to  estimate  the  thermal  diffusivity  of  the  ice.  To  begin,  a  vertical  coordinate 
(z)  is  established  with  origin  at  the  ice/water  interface  and  positive  direction  upward  toward  the 
ice/air  interface.  With  this  coordinate  system,  a  known  ice  thickness  of  42  cm,  and  the  thermistor 
readings  at  the  ice/air  and  ice/water  interfaces,  the  following  initial  and  boundary  conditions  were 
established: 

7'(2,0)  =  -  1.8^C-  (0.5  -)2, 

cm 

7(0./)  =  -l.S-C.  (^2) 


and 


7  (42cm, /)  =  -22.8=C+ (1.54  X  )()■'' - )i.  (A3) 

second 

Equation  (A  1 )  is  a  linear  approximation  of  the  temperature  profile  at  the  beginning  of  the  warming 
trend.  Equation  (A2)  is  the  approximate  water/ice  interface  temperature.  Equation  (A3)  is  a  linear 
approximation  to  the  time  dependence  of  the  upper  ihennisior  during  the  warming  trend. 

After  rewriting  these  conditions  using  the  variable  V  =  /  -  18,  we  solve  the  homogeneous 
equation  of  heat  conduction 

_  !  (7/ 

;,,2  ■  K-rl,  (A41 
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with  the  conditions 


7(2,0)  =  -0.52 


(A5) 


7(0,/)  =  0 


(A6) 


and 


7(42,/)  =  -21.0+  (1.54x10“^)/ 


(A7) 


In  Eq.  (A4)  the  parameter  k  is  the  thermal  diffusivity  who.se  numerical  value  we  shall  predict. 
Next,  Eq.  (A4)  is  Laplace  transformed  giving  (Ref.  21,  p.  420) 


dz 


( 2,  j)  -  - 7 (2,  .0  =  -  ( -  )  7 ( 2,  /  s  0) 
2  K  K 


(A8) 


with  the  transformed  boundary  conditions 


7(0,  ,s)  =  0 


(A9) 


and 


7(42.5) 


21.0  0.54  X  10  ■^) 


(AlO) 


Now  the  object  is  to  determine  7  (z.  5)  and  inverse  iransfomi  to  obtain  /c,  ;t  .  We  proceed 
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by  finding  the  Green  function  which  solves 


=  -5(2-  2') 


(All) 


with  boundary  conditions  g^(O.z')  =  0  •  =  0  .  Standard  techniques  (Ref.  21, 

p.  247)  give 


sinh  f ' 


1/2 


C^C.z')  = 


1/2 


[(^)  zjsmh[(f)  (2^-42) 

^  W  ^1/2 

(-)  sinh(-(-)  42) 

k:  k 


(A  12) 


where  is  the  smaller  of  (2,  2')  and  2^  is  the  larger  of  (2,  z). 

Using  the  Green  function  defined  in  Eq.  (A12),  the  solution  for  7(r,  s)  is 
42 

T(z,s)  =  (—){ z'G  (2',2)d:' -  [G  z) -y  j s)  ~  f  [z',  s)  ~,G  (z'.  z)  “ . 
K  J  ^  dz  dz  ^ 


(A  13) 


0 


0 


Substituting  Eqs.  (A9),  (A  10),  and  (A  12)  into  Eq.  (A  1 3),  one  finds  that 


7'(z..0 


-0.52  fl.54xio' 


I 


-4  , 


) 


f  ^  s  ) 

'  -sinh  I  (  -  )  2  I 

V  K  j 

/  ^  ‘ 

sinh  (-)  42 

V  K'  J 


(A14) 


Inverse  transforming’^  Eq.  (A  14)  and  using  the  definition  of  7  (2,  0  in  terms  of  7(2,0  gives 


,  _  “C  _4 

7  iz,  I )  =  -  1 .8  C  -  ( 0,5  -  )  2  -r  ( 1 .54  X  to  ) 

cm  set  on  it 


IH 


f42(n.  ,) 


42  (n  +  ;,  ) 


>,  [Ku]  '  '  2  (  VC/. )  ’  ■  (Kij)  ‘  ■ 


du , 


-oc() 


(A1.5l 
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FIG.  Al.  Data  points  for  the  last  temperature  profile  taken  during  the  ice  growth  phase  (run  30, 
Fig.  8  )  are  compared  with  theoretical  predictions  based  on  Eq.  (A  15).  Predictions  are  shown  using 
three  different  thermal  diffusivities.  The  best  visual  fit  is  for  a  k  of  0.008  cm^/s. 
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where  erf  'xs,  the  generalized  error  function. 

Equation  (A  15)  was  evaluated  numerically  at  a  time  corresponding  to  the  last  temperature  pro¬ 
file  taken  (run  30,  Fig.  8)  during  the  ice  growth  phase  of  the  experiment.  For  the  parameters  here, 
including  the  terms  (-3)<n<3  was  adequate  for  convergence.  Figure  A1  shows  the  experimental 
temperature  profile  and  the  result  of  evaluating  Eq.  (A  15)  for  three  values  of  thermal  diffusivity  k. 
The  value  of  0.0080  cm^/s  gives  the  best  visual  fit  for  temperatures  less  than  -6°C. 

The  thermal  diffusivity  of  the  ice  is  in  reality  a  function  of  S  (salinity  in  ppt),  T  (temperature 
in  degrees  Celsius),  and  P  (the  density  in  grams  per  cubic  centimeter).  The  equations  in  Ref.  12 
give 


(0.005 


K  = 


(cm)  5(  C)  (cm)5  T 


0.505-^  +  ro.0018— ^ 


s(“C) 


gccy 


src)-i2' 


for  temperatures  above  -8®C.  In  Eq.  (A  16),  k  is  the  thermal  conductivity  of  the  ice  and  c,  is  its  spe¬ 
cific  heat.  With  a  salinity  of  about  7  ppt,  a  density  of  0.91  g/cm^  and  a  temperature  of  -8°C,  Eq. 
(A  16)  predicts  a  diffusivity  of  0.0054  cm”/s.  In  contrast,  pure  ice  has  a  thermal  diffusivity  of 
about  0.012  cm  /s.  The  value  obtained  experimentally  here  is  in  keeping  with  the  limits  estab¬ 
lished  by  these  previous  results. 

After  Submergence 

Here  we  will  use  the  warming  of  the  ice  block  during  the  initial  16  hours  of  submergence  to 
obtain  an  estimate  of  the  thermal  diffusivity.  Equation  (A4)  is  again  applicable,  but  this  time  a  so¬ 
lution  of  the  form 


7'(r,  0  =  R  {r)S{ip)  Z{2)c 


fA17) 
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is  assumed.  Equation  (A  17)  reflects  the  cylindrical  symmetry  of  the  submerged  block;  z  is  now  de¬ 
fined  as  a  down  pointing  vertical  axis  with  2=0  at  the  top  of  the  block,  r  is  the  radial  coordinate 
with  r  =  0  at  the  center  of  the  block,  and  cp  is  the  azimuthal  angle.  The  boundary  conditions  are 


2(0)  =  Z(2„)  =  R(a)  =  0 


(A  18) 


where  zq  is  the  block  thickness  and  a  is  its  radius.  The  angular  symmetry  of  the  boundary  condi- 

1 

tions  imply  no  cp  dependence  in  the  solution.  Using  Eq.  (A18),  a  general  solution  toEq.  (A17)  is“ 


(C.,.)  f  ](  2 


(  — )  sin  (/57t  -  )  e  ,  (A19) 


where  n  and  p  are  summation  indices  going  from  1  to  infinity  and 


p  :  V  a  y  U’n 


(A20) 


In  Eq.  (A  19)  they’s  are  Bessel  functions  and  the  parameters  satisfy  the  relation 


,/o(a  )  =  0 


(A21) 


In  Eq.  (A20)  k  is  the  thermal  diffusivity. 

To  determine  the  constants  ^  in  Eq.  (A  19),  the  initial  conditions  must  be  specified.  Using 
the  experimental  data,  we  have  =  44cm,  n  =  29cm  and,  shortly  after  submergence  (the  first  ther¬ 
mistor  data  was  taken  about  25  minutes  after  submergence),  the  approximate  initial  condition 


7'(r,  z,0)  (-8.7)  (  ,^)//(z)//(  11  -  z)  -  8.7^1  -  1// ( z  -  1  1  )//(:,.-  1  1  >  ' 

11  1,  tz.,  -  1  1  )  ,  "  u 


(A22) 
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In  Eq.  (A22)  the  H's  are  Heaviside  functions  and  it  is  understood  that  at  z  -  1 1  cm  only  one  term 
in  the  bracket  applies.  The  z  dependence  of  the  initial  condition  was  set  by  the  vertical  array  ther¬ 
mistor  data  and  assuming  the  array  was  located  at  the  center  of  the  block  (it  was  actually  a  few  cen¬ 
timeters  to  one  side).  Since  we  don’t  have  a  horizontal  array  of  thermistors,  the  r  dependence  was 
set  such  that  it  has  a  linear  drop  from  the  center  to  the  edge  of  the  block.  This  r  dependence  is  sug¬ 
gested  both  by  the  initial  z  dependence  and  the  symmetry  of  the  forcing  function  (the  water  si.ir- 
rounding  the  block). 

Using  the  initial  condition  in  Eq.  (A22),  the  constants  become 


+  8.7  (; 


n  !  \pn) 


11  ,  f '0 

)  —  ;  cos (pn)  -  cos 


8.7 


(zy  -  1 1)  [pn 

,2 


0 

pn) 

11 {pn) 
^0 


(zn-  11) 


(Zq) 

-  ros{pn) 

pn 


ypn)  [  Zo  J  [pn)  [  z„  I  I 


!  -^1  1  _  r 

i  (a  )  ) 

t_  n  \  I  n  {) 


x^J^^{x)  dx 


(A23j 


Together,  Eqs.  (A19)  and  (A23),  and  the  definition  of  T  in  terms  of  f ,  give  the'solution  for 
the  spatial  and  time  dependence  of  the  block’s  temperature.  These  equations  were  implemented 
and  the  experimental  temperature  profiles  at  f  =  0,  7.6  h,  and  13.6  h  used  to  determine  an  effective 
thermal  diffusivity.  Figure  A2  shows  the  experimental  data  and  calculat’ons  using  two  different 
diffusivities.  The  results  indicate  an  effective  diffusivity  between  0.0013  and  0.0018  cm7s.  .A 
qualitative  explanation  for  this  much  lower  value  of  thermal  diffusivity  (compared  with  the  result 
obtained  in  the  first  part  of  the  appendix)  is  given  in  the  main  body  of  the  paper. 
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wave  speeds  throughout  the  ice  (especially  the  region  near  the  ice/water  interface);  high  frequency 
acoustic  backscatter  near  normal  incidence;  and  their  governing  physical  parameters. 

APPROACH 

The  skeletal  transition  zone  is  a  primary  factor  in  high  frequency  acoustic  characteristics  of  the  ice. 
Under  the  sponsorship  of  ONT,  we  designed  an  experiment  to  measure  the  sound  speed  and  temperature 
profiles  within  growing  and  submerged  sea  ice  while  minimizing  the  disturbance  to  the  skeletal  transition 
zone.  With  this  grant,  we  received  the  funding  necessary  to  analyze  and  publish  the  data  and  its  analysis 
as  wen  as  the  underlying  theory. 

The  experiment  and  the  data  collected  during  March  1990  were  designed  to  allow  examination  of 
vertical  sound  speed  and  temperature  profiles  throughout  the  ice  at  a  resolution  not  previously  obtained 
experimentally.  The  experiment  also  allowed  us  to  relate  these  profiles  with  both  the  salinity  profile 
within  the  ice  as  well  as  with  acoustic  backscaner  from  the  ice. 

In  the  experiment,  linear  arrays  (also  called  strings)  of  ihermisters  and  transducers  were  suspended 
vertically  in  two  large  holes  cut  through  the  ice  canopy  down  to  the  underlying  water.  We  then  allowed 
the  water  to  freeze,  capturing  the  instruments  within  the  growing  ice.  We  then  made  in-situ 
measurements  of  the  acoustic  and  thermal  vertical  profiles  within  the  ice.  (For  the  former,  a  high 
frequency  source  below  the  growing  ice  provided  the  propagating  wave-form  received  at  adjacent 
transducers  to  allow  calculation  of  the  sound  speed  between  transducer  pairs.) 

Immediately  after  each  data  acquisition  session  using  the  receiver  string,  the  transmitter  was  placed 
into  a  transmit/receive  mode  to  measure  normal  incident  backscatter  from  the  ice  at  several  frequencies. 
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We  extracted  effective  interface  reflection  coefficients  from  these  data,  calculated  a  theoretical  sound 
speed  profile  within  the  skeletal  layer  to  compare  with  tiie  measured  profile. 

After  the  ice  had  refrozen  to  where  the  ice/water  interface  was  between  the  deepest  two  receivers 
on  the  linear  arrays,  a  cylindrical  block  was  cut  out  of  the  ice  such  that  the  receiver  and  thermister  strings 
were  centered  in  the  block.  The  block  was  submerged  and,  during  the  hours  it  took  to  come  to 
equilibrium  with  the  water  temperature,  data  similar  to  that  discussed  above  were  acquired.  This  allowed 
monitoring  of  the  changes  in  velocities  and  scattering  properties  in  a  block  under  similar  conditions  to 
blocks  submerged  due  to  ice  rafting  and  ridging.  We  analyzed  these  data  in  the  same  way  as  for  growing 
ice. 

ACCOMPLISHMENTS 

•  Reduction  of  the  data  for  compressional  sound  speed  profiles  and  temperature  profiles  during  the 
ice  growth  and  submergence  stages.  These  have  high  resolution  in  both  time  and  in  the  vertical. 

•  Inferred  bounds  on  and  a  profile  of  attenuation  within  the  ice.  The  skeletal  zone  attenuates  sound 
more  efficiently  than  the  bulk  of  the  ice. 

•  Successful  comparison  of  averages  of  the  directly  measured,  high  resolution  compressional  sound 
speed  profiles  and  that  inferred  from  measured  average  reflection  coefficients  via  our  theory,  for 
growing  ice  (Figure  1)  and  submerged  ice  (Figure  2). 

•  Inferred  thermal  diffusivity  for  growing  and  for  submerged  sea  ice.  The  former  is  consistent  with 
bounds  determined  by  measurements  of  diflfusivity  of  fresh  water  and  the  calculations  of  G. 

Maykut.  The  latter  suggests  that  the  skeletal  layer  that  forms  across  the  entire  surface  of 
submerged  ice  acts  as  an  insulator,  decreasing  the  effective  thermal  diffusivity  by  reducing  heat 
flow  across  the  ice/water  interface. 

•  Inferred  qualitative  brine-drainage  channel  structure  within  the  bottom  one  third  of  growing  ice 
based  upon  the  response  of  the  ice  to  submergence. 

•  Completing  and  will  be  submitting  a  paper  on  the  experiment  and  data  analysis  listed  above. 

•  Submitted  paper  on  theory  for  high  frequency  acoustic  backscatter  from  the  rough  ice/water 
interface  of  first  year  arctic  sea  ice.  Results  include  predictions  of:  sound  speed  profiles  within  the 
skeletal  layer  and  reflection  coefficients  consistent  with  the  measurements  noted  above;  backscatter 
from  ice  blocks,  consistent  with  observations  (Figure  3);  and  one  dimensional  (isotropic)  rough- 
surface  wave-number  spectrum. 

SUMMARY 

We  have  analyzed  the  high-frequency  acoustic  properties  of  the  ice/water  interface  via  direct 
measurement  of  the  in  situ  vertical,  temporal  variation  of  profiles  of  compressional  sound  speed, 
attenuation  and  temperature  of  the  ice  adjacent  to  this  interface.  We  find  that  there  strong  gradients  in  the 
acoustic  properties  of  this  interface  associated  with  gradients  in  physical  properties  of  the  ice.  We  have 
examined  both  the  characteristics  and  reasons  for  these  gradients.  They  affect  the  effective  reflectivity 
and  the  acoustic  backscatter  from  the  arctic  sea  ice  canopy.  We  have  demonstrated  a  theory  for 
backscatter  and  reflectivity  that  successfully  reproduces  the  observations  by  taking  into  account  these 
gradients  within  the  ice  and  the  roughness  of  the  ice  interface. 
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Figure  3.  (A)  IndividuaJ  observed  realizations  of  backscatter  from  an  ice  block  compared  with  theoretical  predic¬ 
tions  of  average  backscatter;  (B)  predicted  compressional  sound  speed  profile  and  (C)  associated  frequency- 
dependent  reflection  coefficient  for  the  ice  used  in  (A). 
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